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A MULTIWAVELENGTH STUDY OF THREE HYBRID 
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ABSTRACT 

We present multiwavelength imaging observations of PKS 1045—188, 8C 
1849+670, and PKS 2216—038, three radio-loud active galactic nuclei from the 
MOJAVE-Chandra Sample that straddle the Fanaroff-Riley (FR) boundary be¬ 
tween low- and high-power jets. These hybrid sources provide an excellent op¬ 
portunity to study jet emission mechanisms and the influence of the external 
environment. We used archival VLA observations, and new Hubble and Chan¬ 
dra observations to identify and study the spectral properties of five knots in PKS 
1045—188, two knots in 8C 1849+670, and three knots in PKS 2216—038. For 
the seven X-ray visible knots, we constructed and fit the broadband spectra us¬ 
ing synchrotron and inverse Compton/cosmic microwave background (IC/CMB) 
emission models. In all cases, we found that the lack of detected optical emission 
ruled out the X-ray emission from the same electron population that produces 
radio emission. All three sources have high total extended radio power, similar 
to that of FR II sources. We find this is in good agreement with previously stud¬ 
ied hybrid sources, where high-power hybrid sources emit X-rays via IC/CMB 
and the low-power hybrid sources emit X-rays via synchrotron emission. This 
supports the idea that it is total radio power rather than FR morphology that 
determines the X-ray emission mechanism. We found no significant asymme¬ 
tries in the diffuse X-ray emission surrounding the host galaxies. Sources PKS 
1045—188 and 8C 1849+670 show significant differences in their radio and X-ray 
termination points, which may result from the deceleration of highly relativistic 
bulk motion. 
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Subject headings: galaxies: active, galaxies: jets, quasars: general, BL Lacertae 
objects: general 


Introduction 


Radio-loud active galactic nuclei (AGNs) encompass a wide range of objects that are 
classified by properties such as their luminosity, jet morphology, and orientation with respect 
to the observer. AGNs with jets that are closely-aligned to our line of sight are called blazars, 
and their alignment results in apparent superluminal motion, relativistic Doppler boosting. 


rapid core variability, and high polarization (Angel & Stockman 1980 Blandford & Konigl 


1979). Blazars comprise the weak optical emission-lined BL Lacertae objects (BL Lacs) 
and strong emission-lined Flat Spectrum Radio Quasars (FSRQs). The Fanaroff-Riley type 


I and II (FRI, FRII; Fanaroff & Riley 1974) classification subdivides radio-loud AGNs 
based on their luminosity and jet morphology: FR I sources have lower luminosities and 
jets with their brightest radio-lobe emission located closer to the AGN core than to the 
jet termination, and FR II sources have higher luminosities and jets with more-distant or 
terminal hotspot emission. In the standard unification scheme, BL Lacs are associated with 


FR Is and FSRQ with FR IIs (Urry & Padovani 1995). Factors that may contribute to the 


FR I and II dichotomy include black hole mass, spin, and accretion rate; jet composition; 
and the external environment (e.g., Meier 1999 Baum et ah 1995 Wold et al. 2007). 


High-resolution Ghandra X-ray Observatory (GXO) and Hubble Space Telescope (HST) 
observations of FRI and FRII radio galaxy jets have revealed systematic differences in 
their broad-band (i.e., radio-optical-X-ray) spectral energy distributions (SEDs). The jets 
of FRI sources display a concave-downward shape in the log — log v plot, and can be 
interpreted as synchrotron emission, while the jets of FRII sources display convex-shaped 
SEDs which can be interpreted as inverse-Gompton (IG) emission from up scattered cosmic- 
microwave background (GMB) photons (the IG/GMB model), or as synchrotron emission 
from a different electron population than the one producing the radio and optical emission 


(Worrall 2009). However, there are major uncertainties associated with the X-ray emission 
mechanisms. Greating ~100 kpc jets with highly energetic synchrotron-emitting electrons 
requires in situ acceleration along the jet because of their short lifetimes (few tens of years), 
and the IG/GMB model requires large bulk relativistic motions (Fjet ~ 5) out to the terminal 
hot spots, which does not appear to be supported by radio observations ( |Bridle et al. 1994). 

The Monitoring Of Jets in Active galactic nuclei with VLBA Experiments (MOJAVE; 


Lister et al. 2009) program provides an extensive set of radio observations of AGNs in the 


northern sky. The original sample consisted of all 135 known AGN with J2000 declination 
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6 > —20°, galactic latitude \b\ > 2.5° and VLBA 15 GHz flux density of at least 1.5 Jy (or 
2 Jy for sources with 6 > 0°) at any epoch between 1994 and 2004. The MOJAVE-Chandra 


sample (MGS; Hogan et ah 2011, Kharb et ah 2012) is a subset that focuses on AGN that 
were the best candidates for X-ray jet emission. It consists of 27 MOJAVE quasars and 
FRH radio galaxies that have extended kiloparsec (kpc) scale 1.4 GHz flux density of at 
least 100 mJy and radio structure at least 3" in extent. The fraction of hybrid sources is 
significantly higher in MOJAVE (6% — 8%; Kharb et ah 2010) than in the FIRST survey 


(1%; Gawroiiski et ah 2006) which was selected for total emission at 1.4 GHz. 


Hybrid radio morphology sources have an FRI-type plume-like jet on one side of the 
radio core and an FRH-type collimated jet with a hot spot on the other. Their total 
radio powers are intermediate between FRIs and FRHs (e.g., Gopal-Krishna fc Wiita||2000). 


These sources exhibit jets with different radio powers and morphologies on either side of the 
accretion disk—black hole system, so this class of source has the potential to become the 
touchstone for ideas put forth to explain the FR dichotomy. 


Hybrid sources in themselves preclude the presence of different central engines or jet 
compositions in the FRI/FRH jets, and point towards external factors driving the FR divide. 
It has been suggested that hybrid sources may simply be lower power FRH sources located 
in asymmetric environments (e.g.. Miller & Brandt 2009). This could cause deceleration to 
occur in the jet encountering the environmental asymmetry and convert it into an FRI jet 
(Gawrohski et al. 2006). For the handful of hybrid sources that have been observed with 
Ghandra, the X-ray jets/knots in some cases have concave-downward shaped SEDs like FRIs 
(e.g., Birkinshaw et al.]|2002 ) and convex-shaped SEDs like FRHs in others (e.g., [Sambruna 


et al. 2008). It is possible that this is a result of small number statistics and/or shallow 


observations (typical exposures < 30 ks). To this end, we have acquired deep Ghandra, 
HST, and VLA observations of three hybrid blazars, viz., PKS 1045—188, 8G 1849-1-670 and 
PKS 2216—038, from the MOJAVE-Ghandra sample (Table [^. 


In this paper we present the results of these multi-wavelength observations. It is struc¬ 
tured as follows: In section 2 we detail our observations and data reduction. In section 3 
we describe the reduced images. In section 4 we discuss the SED generation and model 
behaviour. In section 5 we discuss our results, and in section 6 we summarize our findings. 
Throughout this paper we have adopted a cosmology with Hq = 71 km s“^ Mpc“^, Da = 
0.73, rim = 0.27. 
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2. Observations and Data Analysis 


We carried out Chandra and HST observations of the three blazars in March, July and 
October of 2013 (Table [^. We describe below our multi-wavelength observations made with 
Chandra, HST and VLA. 

Our observations were not simultaneous, especially in the case of the archival VLA 
data. Although these sources can be classihed as blazars, which typically have rapidly- 
variable cores, we are studying the kpc-scale jet knots which are signihcantly downstream 
from the cores. The kpc-scale jet knots are large in extent, and given their observed radio 
brightness temperature of ~ 10® K likely have light crossing times of at least 5000 years, 
neglecting possible beaming effects. So far, the only known case where a kpc-scale jet knot 
shows short-term variability is the knot HST-1 in M87 ( Kovalev et al.||2007 ). This should be 
considered an exceptional case due to M87’s proximity (z = 0.0044) and HST-l’s proximity 
to its core (~ 1", or ~ 80 pc projected). 


2.1. Chandra X-ray Observations 

We carried out the Chandra observations using the AXAF CCD Imaging Spectrome¬ 
ter (ACIS)—S3 chip (which is back-illuminated for low-energy response) in the very faint 
(VFAINT) timed mode. In order to minimise pileup from the bright blazar cores, we used 
the 1/8 subarray mode (frame time = 0.441 s). We specihed a range of roll angles so that 
the charge transfer trail of the blazar core did not contaminate the jet emission. 

In our X-ray observations, only one counter jet was detected and not every approaching 
jet knot seen in prior radio observations was detected. A summary of jet type and knot 
detections is given in Table 

We extracted X-ray spectra using the Chandra Interactive Analysis of Observations 
(CIAO) package version 4.6 with calibration database (CALDB) version 4.5.9. We repro¬ 
cessed the data using the chandra_repro script and hltered the data to include only photons 
with energies in the range 0.5 — 10 keV to account for Chandra calibration and low-energy 
quantum efficiency contamination. The blazar core was bright enough to create a noticeable 
CCD readout streak. Techniques exist to replace the streak with a background spectrum, 
but this was not done to avoid biasing it with the background from either side of the source. 
We used the radio images to determine the sizes and locations of regions for the extraction. 
For the jet extraction, we used background regions on the same side of the source as the 
jet but shifted off axis. For the core extraction, we did not use background regions to avoid 
biasing it with any possible background asymmetry. We then generated source and back- 
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Table 1. Source Properties 


Name 

(1) 

RA (J2000) 

(2) 

Dec (J2000) 

(3) 

Redshift 

(4) 

Scale (kpc/") 

(5) 

Papp (c) 
(6) 

on 

(7) 

PKS 1045-188 

10^48""6L621 

-19°9'35.727" 

0.595 

6.653 

10.51 

10.87 

8C 1849-F670 

18'*49”"16N072 

67°5'41.680" 

0.657 

6.954 

23.08 

4.96 

PKS 2216-038 

22'*18”^52N038 

-3°35'36.879" 

0.901 

7.812 

6.73 

16.9 


Note. — Columns are as follows: (1) Source name (2) Right ascension (J2000) (3) Declination 
(J2000) (4) Distance scale (5) Maximum apparent parsec-scale jet speed in units of the speed of 
light (Lister et al.|2013) (6) Derived maximum parsec-scale viewing angle using 0 = 2 arctan(l//3app) 


Table 2. HST and Chandra Observations 


Source Name 

Telescope 

Obs. Date Exp. Time (s) 

PKS 1045-188 

Hubble 

2013 May 16 

2700 

PKS 1045-188 

Chandra 

2013 Mar 18 

56000 

8C 1849-^670 

Hubble 

2013 Jul 19 

3000 

8C 1849-^670 

Chandra 

2013 Sep 02 

34000 

8C 1849-^670 

Chandra 

2013 Sep 05 

46000 

PKS 2216-038 

Hubble 

2013 Jul 30 

2700 

PKS 2216-038 

Chandra 

2013 Jul 30 

56000 











Table 3. 


Jet Knot Specific Flux Densities 


Name 

Component 

Type 

1.4 GHz 
(mJy) 

4.8 GHz 
(mJy) 

8.4 GHz 
(mJy) 

1537 nm 

(+Jy) 

477 nm 

(^Jy) 

1 keV 
(nJy) 

PKS 1045-188 

Knot A 

FR-II 

26.7 

8.9 

4.2 

< 1.9 

< 0.28 

1.0 


Knot B 

FR-II 

33.6 

12.0 

3.0 

< 9.0 

< 0.55 

2.0 


Knot C 

FR-II 

64.7 

23.9 

13.6 

0.137 

0.152 

0.9 


Knot D 

FR-II 

21.2 

6.6 

< 3.3 

< 0.204 

< 0.136 

< 0.1 


Knot E 

FR-II 

25.0 

6.1 

< 3.3 

0.0825 

0.284 

< 0.1 


Counter Jet 

FR-I 

190.5 

68.3 

< 3.3 

< 0.0087 

< 0.0017 

1.1 

8C 1849+670 

Knot A 

FR-I 

6.1 

2.5 

0.6 

< 0.0054 

< 0.0017 

2.2 


Knot B 

FR-I 

5.6 

3.3 

2.2 

< 0.0054 

< 0.0017 

< 0.2 


Counter Jet 

FR-II 

49.1 

17.1 

12.4 

< 0.0054 

< 0.0017 

< 0.3 

PKS 2216-038 

Knot A 

FR-I 

44.8 

17.4 

9.6 

< 0.0063 

< 0.0016 

0.9 


Knot B 

FR-I 

116.9 

51.3 

47.6 

< 0.0063 

< 0.0016 

1.3 


Knot C 

FR-I 

36.1 

20.5 

2.6 

< 0.0063 

< 0.0016 

0.5 


Note. — In all cases, the knots are in the approaching jet. The optical upper limits for knots A, B, 
and D of PKS 1045—188 suffer from significant external contamination (Section |2.2[). 
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ground spectra and the associated Response Matrix Files (RMFs) and Ancillary Response 
Files (ARFs) using the specextract script. For all jet components, we generated both 
unbinned spectra and spectra with one count per bin. 


We modelled the spectra with an absorbed power law using the XSPEC package. Due to 
the low number of counts in individual jet knots (i.e. tens), including the hydrogen column 
density parameter was tested and did not have a signihcant impact on the £t, so it was 


frozen at values from the HEASARC calculator. We used a C statistic ( 

Cash||1979 

) modihed 

for background subtraction called the W statistic ( 

Wachter et al. 

1979 

), which requires the 


data have at least one count per bin. As verihcation, we also £t the spectra using the ISIS 
package with a maximum likelihood statistic and subplex htting method. The ISIS results 
for both the binned and unbinned data agreed with the XSPEC results for the binned data. 
The hnal £t parameters are listed in Table 


In order to study the galactic environment, we looked for asymmetries in the diffuse 
X-ray emission around each source. We used SAOImage DS9, excluded regions of known 
radio jet emission, and compared the background X-ray counts on each side of the source 
using half-circular regions centered on each jet. Based on the radio emission, we consid¬ 
ered emission within a radius of 23" for PKS 1045—188 and 8C 1849-1-670 and 33" for PKS 
2216—038. We did this comparison on unbinned, unsmoothed 500-10000 eV data hies and 
on data hies that had been hltered to energies of 500-1000, 1000-2000, and 2000-10000 eV. 
We sampled the background throughout the rest of the image using regions of the same 
size to estimate the mean and standard deviation for such regions. In all cases, the back¬ 
ground around the source was elevated above the mean as one would expect near an AGN, 
but none of the diherences between the half-circular regions were greater than 1.5 times 
the standard deviation. The percentage diherences (i.e. for background regions A and B, 
{counts A — counts b) I counts b) that would have been required to be considered signihcant 
(greater than 3a) are 18%, 33%, and 23% for PKS 1045—188, 8C 1849-1-670, and PKS 
2216—038 respectively. Attempts have been made to look at nearby galaxy number densi¬ 
ties and asymmetric galactic interactions (Kharb et ah 2014); however, results have been 
inconclusive due to lack of redshift information to conhrm group or cluster measurement. 


2.2. HST Optical Observations 

We carried out the HST observations using the Wide Field Camera 3 (WFC3) through 
the wideband F160W and F475W hlters (1537 nm and 477 nm central wavelengths, respec¬ 
tively). We loaded pipeline drizzled images into SAOImage DS9 and extracted background- 
subtracted counts using the Chandra-Ed Archive Server analysis commands from the DS9 Vir- 















Table 4. X-ray Spectral Analysis 


Name 

Component 

Nh (lO^'^ cm 

r 

uFy (10 erg s ^ cm 

(1) 

(2) 

(3) 

(4) 

(5) 

PKS 1045-188 

Core 

4.1 

1 jyQ+0.02 

^^-0.03 

457.02;jO;?o 


A 

4.1 

-‘-•'^-0.31 



B 

4.1 

1 Qc:+0-26 

i.OO_o 24 

4-89;1:2i 


C 

4.1 

1 Q-| +0.39 

Q i)o+0.84 
^•^^-0.67 


Counter Jet 

4.1 

1 1 r+0.22 

-'-•■‘-^-0.17 

2 74+0-91 

2-‘+-0.74 

8C 1849-F670 

Core 

5.8 

-'-•^'J-0.02 

773.2i;j|f^ 


A 

5.8 

1 qq+0.15 
-*-•30-0.14 

5.29;j;“ 

PKS 2216-038 

Core 

5.8 

-*-•00-0.03 

422 7Q+11-®’*" 

+^^•'0-11.56 


A 

5.8 

1 9C+0.27 

-*-•20-0.24 

9 9/1+0.96 

2-^+-0.74 


B 

5.8 

1 44+0-26 
-*-•++-0.23 

O-1+1-0.89 


C 

5.8 

1 -] ^+0.36 
-'-•-‘-‘-0.31 

1 14+0.71 
-1--1+-0.49 


Note. — Columns are as follows; (1) Source name (2) Region of spectrum extraction (3) 
Hydrogen column density (4) Photon index of the power law distribution (5) at 1 keV. 
The hydrogen column density was fixed at galactic values, and the 90% confidence values are 
given for the other parameters. 
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tual Observatory. We then computed XFx values by multiplying the background-subtracted 
counts (in units of electron s“^) by the inverse sensitivity keyword PHOTFLAM (ergs cm“^ 
electron”^) and the pivot wavelength keyword PHOTPLAM (A) from the image headers. 

For the sources 8C 1849-1-670 and PKS 2216—038, we detected no significant optical 
emission other than the jet core. This lack of detection is meaningful for constraining the 
X-ray emission mechanism, so three times the standard deviation of the background was 
used as an upper limit for the SED data points. 

For PKS 1045—188 jet knots A, B, and D, any possible optical emission was contami¬ 
nated by multiple stars and galaxies in the held of view. Because a perfect subtraction was 
impossible for such a crowded held and the optical huxes were vital for ruling out potential 
emission models, we took a conservative approach and used direct, unmodihed measurements 
as upper limits for the SED data points. For PKS 1045—188 knot C, we used a smaller region 
of 0.6" radius to reduce contamination because beyond this radius, when looking in directions 
away from contamination (i.e., away from the galaxy that overlaps the jet between knots B 
and C), count rates fell back to background levels. 


PKS 1045—188 knots C and E have possible optical detections. Small regions of emission 
coincide well with the peaks of the 4.8 GHz radio knots. To estimate the likelihood of this 
being coincidence, we performed a simple Poisson probability test. We used the daof ind task 
in the DAOPHOT ( Stetson||1987 ) package in IRAF (Image Reduction and Analysis Facility) 
version 2.16 to count the number of stars in the held, and then used Poisson statistics 
to estimate the likelihood of hnding a random source with equal or lesser magnitude in the 
knot regions. We used regions of 0.5" radius based on the Chandra PSF size. The chance 
probabilities are 4.4% and 2.1% for knot C and E respectively, which supports these being 
jet detections. 


2.3. VLA Radio Observations 

Pipeline calibrated radio data for all three blazars were available at several frequencies 
and array conhgurations of the VLA in the NRAO image archiv^ We specihcally examined 
the VLA A-array data at 1.4, 4.8 and 8.4 GHz (Table [^. We downloaded the calibrated 
{u,v) datasets and after additional phase self-calibration, created hnal images convolved with 
the appropriate beam-sizes to match the Chandra and HST data. Imaging and further self- 


%ttp://iraf.noao.edu/ 

^https://archive.nrao.edu/archive/archiveimage.html 
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calibration were carried out using the Difmap package (Shepherd 1997). We began with a 
human-guided CLEAN deconvolution to ensure that CLEAN windows coincided with regions of 
actual emission, and once the jet knots were identified, further iterations of CLEAN-ing and 
self-calibration were done using a loop until the model converged. To estimate the fluxes of 
individual features, we replaced the CLEAN components in the region of the jet knots with 
circular Gaussian models using the Difmap modelfit task and then used the model error 
analysis tool in the Difwrap package (Lovell 2000). In the event of no detection, we took 
three times the standard deviation of the background as an upper limit. 


3. Image Analysis 
3.1. PKS 1045-188 


Given that MOJAVE VLBA imaging shows a one-sided parsec-scale jet at position angle 
146° (measured from north through east), we consider the kiloparsec-scale jet at position 
angle 125° (Lister et ah 2013) to correspond to the approaching jet. The 4.8 GHz VLA 
image shows five FRII approaching jet knots and the FRI counter jet plume. Extensive 
X-ray emission is present along the approaching jet up to and including knot G, after which 
the radio jet bends ~90° and continues to knots D and E (Figure [^. The radio jet bend is 
likely exaggerated by the small angle between the jet and line of sight {6), which could well 
be comparable to the angle of the parsec-scale jet to the line of sight, 9 <2 aictaxi{l/(3max) = 
10.9°. There is possible X-ray emission coincident with the radio counter jet plume, but it is 
too faint and diffuse to concretely identify it as counter jet emission without deeper exposure. 
Any possible optical emission from knots A and B is obscured by or blended with two field 
galaxies. There is faint optical emission at the sites of knots G and E, but no optical emission 
is present at the locations of knot D or the counter jet lobe (Figures]^ and |^. 


3.2. 8C 1849+670 


Given that MOJAVE VLBA imaging shows a one-sided parsec-scale jet at position 
angle —52°, we consider the kiloparsec-scale jet at position angle 0° (Lister et ah 2013) to 


correspond to the approaching jet. The 1.4 GHz VLA image shows two FRI approaching jet 
knots and the terminal FRII counter jet hotspot. X-ray emission is detectable only from the 
approaching jet knot A (Figure]^. No optical emission is detected from either jet (Figures]^ 
and|^. 
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3.3. PKS 2216-038 


Given that MOJAVE VLBA imaging shows a one-sided parsec-scale jet at position 
angle —170°, we consider the kiloparsec-scale jet at position angle 135° ( [Lister et ah 2013) 
to correspond to the approaching jet. The 1.4 GHz VLA image shows three FRI approaching 
jet knots and the FRII counter jet hotspot. Extensive X-ray emission is present throughout 
the approaching jet (Figure]^. The HST image shows a field galaxy near the optical core 
position, but no clear jet emission (Figures]^ and |^. It is not clear if this galaxy is physically 
close to the host galaxy of PKS 2216—038 because no spectroscopic information is available 
on it. 


4. SED Generation 


We used the jet emission model of Krawczynski et ah (2004) to construct fits to the 
jet knot SEDs. It generates spectral energy distributions for synchrotron, inverse-Gompton 
(IG) of cosmic microwave background (GMB) photons, and synchrotron self-Gompton (SSG) 
emission for a given set of physical jet parameters. We used magnetic field and bulk flow 
Doppler factor values calculated from VLBI observations by Hogan et al.| (2011 ) as a starting 
point and adjusted them to fit. The SEDs are presented in Figures [T0| [TTt [T^and the SED 
best-fit parameters are given in Table 


The jet parameters include Doppler factor (J), magnetic field strength (B), electron 
energy density minimum and maximum electron Lorentz factors {'jmin and 'fmax), 

and power-law index of the electron energy distribution (n). The Doppler factor, electron 
energy density, and minimum Lorentz factor collectively affect the predicted fluxes of all three 
spectral components, but the maximum Lorentz factor has little influence on the predicted 
fluxes. The magnetic field also affects the predicted fluxes of the synchrotron and SSG 
components. The minimum and maximum Lorentz factors determine the range of frequency 
values spanned by all three spectral components. 


We constructed best-fit SEDs for all knots detected in X-rays, which resulted in three 
SEDs for PKS 1045—188, one for 8G 1849-1-670, and three for PKS 2216—038. Gore SEDs 
were not created because of optical saturation and contamination. In all cases, the SSG 
component is several orders of magnitude too low to fit the knot fluxes, but the synchrotron 
and inverse-Gompton components provide a reasonable fit for the radio and x-ray compo¬ 
nents, respectively. The best fit model parameters are not well constrained, as different 
combinations of parameters could achieve the same fits. With additional observations, the 
spectral index could be precisely determined, and with observations at the highest and lowest 
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Table 5. VLA Observations 


Name 

Frequency 

(GHz) 

Project 

ID 

Obs. 

Date 

Beam 
(" X ") 

Position Angle 

n 

Image rms 
(mJy beam“^) 

PKS 1045-188 

1.49 

AG361 

1992 Nov 18 

3.18x1.29 

-29.36 

0.27 


4.86 

AB660 

1992 Dec 14 

1.13x0.38 

-42.38 

0.33 


8.44 

AG361 

1992 Nov 18 

0.47x0.16 

10.27 

0.17 

8C 1849-F670 

1.51 

AL499 

1999 Aug 27 

1.43x1.18 

-8.53 

0.17 


4.86 

AP250 

1993 May 02 

1.54x1.38 

-62.17 

0.09 


8.44 

AL401 

1997 Aug 26 

3.19x2.42 

4.42 

0.11 

PKS 2216-038 

1.40 

AL634 

2004 Nov 21 

0.41x0.54 

27.01 

0.17 


4.86 

AK491 

1999 Aug 05 

5.02x3.92 

17.86 

0.17 


8.46 

AR415 

1999 Jul 31 

3.82x2.36 

39.67 

0.62 



07.00s 10h48m06.00s 

RA {J2000) 



07.00s 10h48m06.00s 

RA (J2000) 


Fig. 1.— (left) Chandra X-ray image of PKS 1045—188 in colonr snperimposed by VLA 4.8 
GHz radio contonrs with contrast adjnsted to show the possible connter jet emission (Section 
3.1[ ). (right) Jet knots identified. The color scales correspond to image connts. The lowest 
contonr level is three times the radio image rms, and each higher contonr is fonr times the 


previons one. 
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RA {J2000) RA (J2000) 

Fig. 2.— VLA 4.8 GHz radio contours superimposed on (left) the colour HST/F160W 
image and (right) HST/F475W image of PKS 1045—188. The color scales correspond to 
image counts. The lowest contour level is three times the radio image rms, and each higher 
contour is four times the previous one. 



RA (J2000) RA (J2000) 

Fig. 3.— Chandra X-ray contours superimposed on (left) the colour HST/F160W image and 
(right) HST/F475W image of PKS 1045—188. The color scales correspond to image counts. 
The lowest contour level is three times the radio image rms, and each higher contour is four 
times the previous one. 
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Fig. 4.— (left) Chandra X-ray image of 8C 1849-1-670 in colour superimposed by VLA 4.8 
GHz radio contours, (right) Jet knots identihed. The color scales correspond to image 
counts. The lowest contour level is three times the radio image rms, and each higher contour 
is four times the previous one. 



RA (J2000) RA (J2000) 

Fig. 5.— VLA 4.8 GHz radio contours superimposed on (left) the colour HST/F160W image 
and (right) HST/F475W image of 8G 1849-1-670. The color scales correspond to image counts. 
The lowest contour level is three times the radio image rms, and each higher contour is four 
times the previous one. 
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RA (J2000) RA (J2000) 


Fig. 6.— Chandra X-ray contours superimposed on (left) the colour HST/F160W image and 
(right) HST/F475W image of 8C 1849-1-670. The color scales correspond to image counts. 
The lowest contour level is three times the radio image rms, and each higher contour is four 
times the previous one. 
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Fig. 7.— (left) Chandra X-ray image of PKS 2216—038 in colour superimposed by VLA 
1.4 GHz radio contours, (right) Jet knots identihed. The color scales correspond to image 
counts. The lowest contour level is three times the radio image rms, and each higher contour 
is four times the previous one. 
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Fig. 8.— VLA 1.4 GHz radio contours superimposed on (left) the colour HST/F160W 
image and (right) HST/F475W image of PKS 2216—038. The color scales correspond to 
image counts. The lowest contour level is three times the radio image rms, and each higher 
contour is four times the previous one. 
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Fig. 9.— Chandra X-ray contours superimposed on (left) the colour HST/F160W image and 
(right) HST/F475W image of PKS 2216—038. The color scales correspond to image counts. 
The lowest contour level is three times the radio image rms, and each higher contour is four 
times the previous one. 
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Fig. 10.— Broad-band spectral energy distribntions (SEDs) of (top) knot A, and (middle) 
knot B, (bottom) knot C of PKS 1045—188. The solid bine line, red dotted line, and green- 
dotted-dashed line are synchrotron, IC/CMB, and SSC spectral components, respectively. 
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Fig. 11.— Broad-band SEDs of knot A in the 8C 1849-1-670 jet. The solid blue line, 
red dotted line, and green-dotted-dashed line are synchrotron, IC/CMB, and SSC spectral 
components, respectively. 


Table 6. SED best-ht parameters 


Name 

Knot 

5 

B 

'^Psoll 

'Tmm 

'^max 

n 




(pG) 

(10“^^ erg cm“^) 


(10^) 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

PKS 1045-188 

A 

4.1 ± 0.7 

45 + 6 

5 + 12 

12 + 17 

1 + 9 

4.0 + 0.1 


B 

4.1 ± 0.7 

50 + 8 

13 + 29 

14 + 19 

2.5 + 9 

4.2 + 0.2 


C 

4.1 ± 0.7 

70 + 10 

4 + 13 

12 + 2 

2.5 + 14 

4.0 + 0.1 

1849+670 

A 

4.2 ± 0.7 

10 + 4 

0.35 + 2 

30 + 14 

12 + 10 

3.5 + 0.2 

PKS 2216-038 

A 

3.0 ± 0.9 

20 + 7 

0.19 + 3 

26 + 16 

3.5 + 2 

3.1 + 0.3 


B 

3.0 ± 0.8 

50 + 3 

0.29 + 3 

32 + 17 

2 + 2 

3.4 + 0.1 


C 

3.0 ± 0.9 

15 + 4 

0.1 + 1 

20 + 19 

3.3 + 2 

2.7 + 0.3 


Note. — Columns are as follows: (1) Source name (2) Region of spectrum extraction and htting 
(3) Bulk flow Doppler factor (4) Magnetic held strength (5) Electron energy density (6) Minimum 
electron Lorentz factor (7) Maximum electron Lorentz factor (8) Power-law index of the electron 
energy distribution. Errors given are one standard deviation from random sampling. 
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Fig. 12.— Broad-band SEDs of knots A, B and C in the PKS 2216—038 jet. The solid 
bine line, red dotted line, and green-dotted-dashed line are synchrotron, IC/CMB, and SSC 
spectral components, respectively. 
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wavelengths of the distributions, the electron Lorentz factors could be better constrained. 
Unfortunately there are very few instruments with the requisite sensitivity and angular res¬ 
olution capability that can obtain such data; therefore, most jet SED studies such as ours 
are currently limited in this way. 

We estimated uncertainties in the £t parameters factors through sampling. We tested 
randomly chosen sets of parameters until we found 100 sets that successfully fit the ob¬ 
served data, and then we calculated the standard deviation for each parameter. Reported 
uncertainties are one standard deviation. 


5. Discussion 


Although the SEDs do not uniquely determine the physical jet parameters, the low 
optical fluxes rule out the pure synchrotron emission for radio through X-ray. For PKS 
1045—188 knots A and B, the optical upper limits are just high enough to allow a pure 
synchrotron model, but we stress that the majority of the fluxes of those upper limits likely 
comes from external contamination (Section |2.2 ). We hnd that the X-ray emission from the 
jets of all the three blazars is consistent with the IC/CMB mechanism, but due to the small 
number of spectral points other X-ray emission models that do not contribute to the optical 
emission may work as well. In all three sources, the X-rays are found on the side with the 
approaching radio jet. The X-rays are detected on the FRI side (i.e., the side without a 
terminal hot spot) in 8C 1849-1-670 and PKS 2216—038, but on the FRII side (the side with 
a radio hot spot, e.g. knot C) in PKS 1045—188. 

Although there are no detectable differences in the diffuse X-ray environment (Section 


2.1), for PKS 1045—188 and 8C 1849-1-670 in both cases the X-ray jets terminate signihcantly 


earlier than their corresponding radio jets (Fig. [^and|^. In addition, for PKS 1045—188, 
the radio jet bends signihcantly at the X-ray terminal point. The X-ray terminal peaks may 
indicate the sites of bulk jet deceleration, where the jet Doppler factor and the IC/CMB 
emission both drop dramatically. 

We searched the literature for hybrid/asymmetric radio sources with Chandra data and 
found discording results as far as the X-ray emission mechanism was concerned (Table [^. In 


only one case, 3C371 (Sambruna et ah 2007), the X-ray emission was detected on the FRII 


side and was consistent with synchrotron emission. In all other cases, the X-ray emission 
was detected on the FRI side. It was consistent with synchrotron emission in 0521—365 
(Birkinshaw et ah 2002), 2201-1-044 (Sambruna et ah 2007[ ), NGC6251 (especially the inner 


jet; Evans et al. 2005), and it was consistent with IC/CMB in PC 1004-1-130 (Miller et ah 
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2006), 3C433 and 4C65.15 (Miller & Brandt 2009), 2007+777 (Sambrnna et al. 2008[ ), and 


3C 17 (Massaro et al. 2009). It is not clear if the X-ray emission is synchrotron or IC/CMB 
in the FRI jet of 1510—089 (Sambrnna et al. 2004). While the X-ray jet in this sonrce is 
mnch shorter (less than half in extent) than the radio jet, similar to synchrotron jets in FRI 
sources, the lack of a clear optical counterpart favours the IC/CMB mechanism. 


In 11 of the 13 studied hybrid sources, the FRI jet is the approaching. X-ray emitting 
jet. This trend could be a coincidence due to the small-number statistics, but if holds true 
after studies of additional hybrid sources then it could seriously question theories regarding 
the creation and the true nature of hybrid sources. Those theories that rely on asymmetries. 



be the approaching jet in the majority of cases. Other theories would remain, such as if 
hybrid sources were actually relativistic FRII sources with bent approaching jets, then the 
inner portion of the jet could be beamed toward the observer and the terminal hotspot could 
be beamed in a different direction, giving the approaching jet the appearance of an ’’edge- 
darkened” FRI. Again, we must emphasize that until statistics are improved, nothing can 
be conclusively ruled out, so studies of additional hybrid sources are vital. 


Our observations indicate that even jets with low power FRI-like appearance have X- 
rays from IC/CMB when the total radio power of the source is like an FRII. This could 
imply that even on the FRI side there is a fast collimated inner spine, just like an FRII 
jet, but which dissipates or is de-beamed before a terminal hot spot is formed. It appears 
that, rather than the FR morphology, the total radio power is the main determinant for the 
X-rays being IC/CMB (in high power sources) or synchrotron (in low power sources). This 
is supported by the plot of total radio power versus redshift (Figure [I^. With the exception 
of 2007+777 (which is a borderline FRI/II source), all sources with FRII-like total radio 
powers at 1.4 GHz have X-ray emission from the IC/CMB mechanism. In the remaining 
sources, synchrotron emission dominates. Furthermore, in all of the hybrid sources studied, 
regardless of emission mechanism, the X-ray emission is from the approaching jet. This is in 
good agreement with previous studies, which indicate that blazars are excellent candidates 
for X-ray jet detections (Hogan et al. 2011). 


6. Conclusions 

In this paper we present multiwavelength jet observations of PKS 1045—188, 8C 1849+670, 
and PKS 2216—038, three radio-loud AGN from the MOJAVE-Chandra Sample that strad¬ 
dle the Fanaroff-Riley boundary. These hybrid sources are of interest because they provide an 
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Fig. 13.— Total radio power at 1.4 GHz plotted against redshift for all hybrid jet AGNs 
in the literatnre. The red symbols denote sonrces with X-ray jets consistent with the syn¬ 
chrotron mechanism, and the black symbols denote those consistent with IG/GMB. The 
open stars denote the three sonrces presented in this paper. The dashed lines demarcate FR 
classihcation by Ledlow & Owen (1996). The magenta line indicates the detection limit of 


the VLSS. 
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excellent opportunity to study the jet emission mechanisms and the influence of the external 
environment. 


We used archival VLA observations at 1.4, 4.6, and 8.4 GHz; deep HST observations 
with the F160W and F475W Alters; and deep Chandra observations. We used MOJAVE 
VLBA parsec-scale observations to identify the approaching jets and VLA kiloparsec-scale 
observations to classify the approaching jets as FRII for PKS 1045—188 and FRI for 8C 
1849-1-670 and PKS 2216—038. From the Chandra observations, we identified X-ray emission 
from three of five radio-visible knots in PKS 1045—188, one of two knots in 8C 1849-1-670, 
and all three knots in PKS 2216—038. We also identified X-ray emission from the area of the 
counter jet of PKS 1045—188, but it is too faint and diffuse to identify any specific emission 
regions without further observation. 


For the seven X-ray visible knots, we constructed and fit SEDs using the synchrotron 
and IC/CMB emission models of Krawczynski et ah (2004). Although we would require 
observations at additional wavelengths to fully constrain the model parameters, we found 
that the weak optical emission ruled out synchrotron emission for radio to X-ray in all cases, 
even for the FRI jets of 8C 1849-1-670 and PKS 2216—038. 


All three sources have high total extended radio power, similar to that of FRII sources. 
We And this is in good agreement with previously studied hybrid sources, where high-power 
hybrid sources emit X-rays via IC/CMB and the low-power hybrid sources emit X-rays 
via synchrotron emission. This supports the idea that it is total radio power rather than 
FR morphology that determines the X-ray emission mechanism. Additionally, the X-ray- 
emitting jet is the approaching jet in all three sources, which is also in good agreement with 
previously studied hybrid sources. 


We found no significant asymmetries in the X-ray environments. Sources PKS 1045—188 
and 8C 1849-1-670 show significant differences in their radio and X-ray termination points 
which may be locations of bulk deceleration. 

Questions remain regarding why the FRI jet is the approaching. X-ray-emitting jet in 
11 of the 13 studied hybrid sources. This could still be coincidence due to the small number 
of studied hybrid sources, or this could indicate additional effects, such as the possibility that 
some of the FRI jets are actually bent, de-beamed FRII jets. Observations of additional 
sources are necessary to answer these questions. 
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(SAO) for and on behalf of NASA under contract NAS8-03060. Support for program number 
13116 was provided by NASA through a grant from the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research in Astronomy, Inc., un¬ 
der NASA contract NAS5-26555. The MOJAVE program is supported under NASA-Fermi 
grant NNX12A087G. The National Radio Astronomy Observatory is a facility of the Na¬ 
tional Science Foundation operated under cooperative agreement by Associated Universities, 
Inc. 

FactUUes: VIA, HST (WFC3), CXO (ACIS). 

REFERENCES 

Angel, J. R. P., & Stockman, H. S. 1980, ARA&A, 18, 321 
Baum, S. A., Zirbel, E. L., & O’Dea, C. P. 1995, ApJ, 451, 88 
Birkinshaw, M., Worrall, D. M., & Hardcastle, M. J. 2002, MNRAS, 335, 142 
Blandford, R. D., & Konigl, A. 1979, ApJ, 232, 34 

Bridle, A. H., Hough, D. H., Lonsdale, C. J., Burns, J. O., & Laing, R. A. 1994, AJ, 108, 
766 

Cash, W. 1979, ApJ, 228, 939 

Evans, D. A., Hardcastle, M. J., Croston, J. H., Worrall, D. M., & Birkinshaw, M. 2005, 
MNRAS, 359, 363 

Fanaroff, B. L., & Riley, J. M. 1974, MNRAS, 167, 31P 

Gawrohski, M. P., Marecki, A., Kunert-Bajraszewska, M., & Kus, A. J. 2006, A&A, 447, 63 
Gopal-Krishna, & Wiita, P. J. 2000, A&A, 363, 507 

Hogan, B. S., Lister, M. L., Kharb, P., Marshall, H. L., & Cooper, N. J. 2011, ApJ, 730, 92 

Kharb, P., Lister, M. L., & Cooper, N. J. 2010, ApJ, 710, 764 

Kharb, P., Lister, M. L., Marshall, H. L., & Hogan, B. S. 2012, ApJ, 748, 81 

Kharb, P., Stanley, E., Lister, M., Marshall, H., O’Dea, C., & Baum, S. 2014, ArXiv e-prints 



- 25 - 

Kovalev, Y. Y., Lister, M. L., Homan, D. C., & Kellermann, K. I. 2007, ApJ, 668, L27 

Krawczynski, H., Hughes, S. B., Horan, D., Aharonian, F., Aller, M. F., Aller, H., Boltwood, 
P., Buckley, J., Coppi, P., Fossati, G., Gotting, N., Holder, J., Horns, D., Kurtanidze, 

O. M., Marscher, A. P., Nikolashvili, M., Remillard, R. A., Sadun, A., & Schroder, 
M. 2004, ApJ, 601, 151 

Ledlow, M. J., & Owen, F. N. 1996, AJ, 112, 9 

Lister, M. L., Aller, M. F., Aller, H. D., Homan, D. G., Kellermann, K. L, Kovalev, Y. Y., 
Pushkarev, A. B., Richards, J. L., Ros, E., & Savolainen, T. 2013, AJ, 146, 120 

Lister, M. L., Gohen, M. H., Homan, D. G., Kadler, M., Kellermann, K. L, Kovalev, Y. Y., 
Ros, E., Savolainen, T., & Zensus, J. A. 2009, AJ, 138, 1874 

Lovell, J. 2000, in Astrophysical Phenomena Revealed by Space VLBI, ed. H. Hirabayashi, 

P. G. Edwards, & D. W. Murphy, 301-304 

Massaro, F., Harris, D. E., Ghiaberge, M., Grandi, P., Macchetto, F. D., Baum, S. A., 
O’Dea, G. P., & Gapetti, A. 2009, ApJ, 696, 980 

Meier, D. L. 1999, ApJ, 522, 753 

Miller, B. R, & Brandt, W. N. 2009, ApJ, 695, 755 

Miller, B. P., Brandt, W. N., Gallagher, S. G., Laor, A., Wills, B. J., Garmire, G. P., & 
Schneider, D. P. 2006, ApJ, 652, 163 

Sambruna, R. M., Donato, D., Gheung, G. G., Tavecchio, F., & Maraschi, L. 2008, ApJ, 684, 
862 

Sambruna, R. M., Donato, D., Tavecchio, F., Maraschi, L., Gheung, G. G., & Urry, G. M. 
2007, ApJ, 670, 74 

Sambruna, R. M., Gambill, J. K., Maraschi, L., Tavecchio, F., Gerutti, R., Gheung, G. G., 
Urry, G. M., & Ghartas, G. 2004, ApJ, 608, 698 

Shepherd, M. G. 1997, in Astronomical Society of the Pacihc Gonference Series, Vol. 125, 
Astronomical Data Analysis Software and Systems VI, ed. G. Hunt & H. Payne, 77 

Stetson, P. B. 1987, PASP, 99, 191 

Urry, G. M., & Padovani, P. 1995, PASP, 107, 803 



Wachter, K., Leach, R., & Kellogg, E. 1979, ApJ, 230, 274 

Wang, J. C. L., Sulkanen, M. E., & Lovelace, R. V. E. 1992, ApJ, 390, 46 

Wold, M., Lacy, M., & Armus, L. 2007, A&A, 470, 531 

Worrall, D. M. 2009, A&A Rev., 17, 1 


This preprint was prepared with the AAS MLX macros v5.2. 



27 


Table 7. Summary of Hybrid Source X-ray Jets 


Name X-ray Jet Type Emission Mechanism Reference 


0521-365 

FRI 

Synchrotron 

2201+044 

FRI 

Synchrotron 

NGC 6251 

FRI 

Synchrotron 

PG 1004+130 

FRI 

IG/GMB 

3G 433 

FRI 

IG/GMB 

4G 65.15 

FRI 

IG/GMB 

2007+777 

FRI 

IG/GMB 

3G 17 

FRI 

IG/GMB 

1510-089 

FRI 

IG/GMB 

8G 1849+670 

FRI 

IG/GMB 

PKS 2216-038 

FRI 

IG/GMB 

3G 371 

FRII 

Synchrotron 

PKS 1045-188 

FRII 

IG/GMB 



This paper 
This paper 


Sambruna et ah (2007) 


This paper 

































































